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Abstract 

Due to e + e~-pair production in the field of supercritical (Z Z cr ~ 
170) nucleus an electron shell, created out of the vacuum, is formed. The 
distribution of the vacuum charge in this shell has been determined for 
super-charged nuclei Ze 3 > 1 within the framework of the Thomas-Fermi 
equation generalized to the relativistic case. For Ze 3 S> 1 the elec- 
tron shell penetrates inside the nucleus and almost completely screens 
its charge. Inside such nucleus the potential takes a constant value equal 
to Vo = — (37r 2 7ip) 1/ ' 3 ~ — 2m-n-c 2 , and super-charged nucleus represents 
an electrically neutral plasma consisting of e,p and n. Near the edge of 
the nucleus a transition layer exists with a width A « a>~ 1 / 2 H/m n c ~ 15 
fm, which is independent of Z (h/mnC <C A <C h/m e c). The electric field 
and surface charge are concentrated in this layer. These results, obtained 
earlier for hypothetical superheavy nuclei with Z ~ A/2 < 10 4 -f- 10 6 , 
are extrapolated to massive nuclear density cores having a mass num- 
ber A » (mpi anc k/m n ) ~ 10 57 . The problem of the gravitational and 
electrodynamical stability of such objects is considered. It is shown that 
for A > 0.04(Z/A) 1//2 (mp iancfc /m n ) 3 the Coulomb repulsion of protons, 
screened by relativistic electrons, can be balanced by gravitational forces. 
The overcritical electric fields E ~ m 2 c 3 /efi are present in the narrow 
transition layer near the core surface. 

The Dirac equation for an electron in the field of a point-like electric charge 
Ze loses its sense for Z > 137, since the energies e n of the bound states ns\/2 
and npi/2 become complex [l]-[3]. For instance, in the case of the lowest energy 
levels one has 



e (l s i/2) = m e c 2 \/l — Q 2 — for the ground state, (1) 



e(2s 1/2 ) = e(2pi /2 ) = m e c 2 



(10 



e(2p 3/2 ) = m e cVl " C 2 /4, -, 
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where < ( = Za < 1, a = e 2 /he = 1/137, m e is the electron mass and the 
potential is assumed to be V(r) = —<^/r, < r < oo. Analogous singularities 
at £ = 1 appear for all ns 1 / 2 and npi/2 states: 

, 2 n-l + ^l-C 

e n /m e c = - 



N 2 +2(n- l)^! 3 ^] 172 

^ + ^3Vl-C + ^5 ) + 

+ 0((W 2 ) 3/2 ), C^l, (2) 



where A = yjn 2 — 2n + 2 and e„ = m e c 2 yl — iV -2 at C = 1- In particular, for 
the case of highly-excited, n 3> 1, states 

e/m e c 2 = l-±- ^ - + < ( < 1. (2') 

' 2n 2 n 3(i+ ^1-C 2 ) 

The appearance of complex energies £ n (C) at ^ > 1 contradicts to unitarity 
and hermiticity of the Dirac Hamiltonian, so an immediate analytical conti- 
tuation of the previous formulae to Z > 137 region is impossible. Analogous 
singularities exist for other physical quantities, for example, for the mean radius 
and the magnetic moment of the ground state, 



(0 = 1 + 2 y^ -*c, ^= J(1 + 2 V / W 2 "W, (1") 

where lc = h/m e c = aag and (1b = eh/2m e c. 

Many aspects of Quantum Electrodynamics of strong fields are considered in 
refs.[5-29], including the relativistic Coulomb problem with Z > 137 [5,6,11-15], 
the critical nuclear charge Z cr [5,11,15], vacuum polarization and superbound 
electrons in the lower continuum at Z > Z cr [10,13,17], spontaneous production 
of positrons at Z > Z cr [6,8,12] and in collisions of two heavy nuclei with 
united charge Z\ + Z 2 > Z cr [10,13,14], the Thomas-Fermi approach for super- 
charged nuclei [7,18], e + e _ -pair creation from vacuum in strong electric field 
and by intense laser radiation [22-27], etc. For further details see the reviews 
[9,10,13,21,28]. 

The difficulty arising at £ > 1, which sometimes is called "Z = 137 catas- 
trophe" [28], is removed if one takes into account the finite size of nucleus [4] 

V(r) = - — f(r/r N ), (3) 
r 



1 Similar situation takes place in other problems with the so-called "fall down to the center" 
in quantum mechanics, see § 35 in ref . [30] . For solutions of the Dirac equation with a point-like 
Coulomb potential "fall down to the center" begins at f = j + 1/2, where j = 1/2, 3/2, 5/2, ... 
is the total angular momentum of the electron state [10,31]. 
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where the cut-off function / depends on the distribution of electric charge inside 
the nucleus. For example, the function 



x(3-x 2 )/2, 0<x<l 
/(,■) = { (3') 
1, x = r/r N > 1 

corresponds to uniform volume density of electric charge Ze and is frequently 
used in calculations. Here the potential at the centre of nucleus is finite, V(0) = 
— 1.5Za/rN. 

When a finite radius rjy > is introduced, the ground level lsi/2 is going 
down monotonously while Z increasing and reaches the boundary of the lower 
continuum e = —m e c 2 for Z = Z cr « 170 [5,11,15]. It can be shown that the 
asymptotic expressions for lsi/ 2 energy are [12,13] 



e(lsuo) { v/W^oth (Av/W^), < C < 1, 

> (4) 

Tc^Tctg (AVc*=T), C>i, 



m e c 2 



where ( — Za, A = ln(^c/rjv) >• 1 is a logarithmic parameter in the problem 
considered and lc — h/m e c = 386 fm is the Compton wave length for electron. 
Eqs.(4) explicitly show that £ — 1 is not a singular point for the energy e((), on 
the contrary to the case of a point-like charge, and energy levels e(n,j) of the 
bound states Is, 2p, 2s, ... smoothly continue to drop into the lower continuum as 
Z increases, until at Z = Z cr (n,j) its energy reaches the boundary of the lower 
continuum. Numerical values of the "critical nuclear charge" Z cr (n,j) were 
obtained by different calculation methods from the equation e(n, j) = —m e c 2 , 
see [5,10,11-13] and references therein. A simple asymptotic formula for Z cr 
follows from eq.(4) and analogous equations for 2p, 2s, ... states: 

where n r = n for ns-states, n r = n — 1 for np^-statcs, n = 1,2,3, ... is the 
principal quantum number, c n j — 2n for ns-states, c n j = \p2 — 1 for 2p!/ 2 -state 
and the condition A 3> n was assumed. As can be seen from the following 
Table, the approximation (5) is rather good for the lowest levels of the electron 
spectrum, though the expansion parameter A <~ 3.7 is not very large. 

Critical nuclear charge for the low-lying states of electron spectrum 
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Atomic 
state 


7(0) 


Z cr 


Z cr 


Ccr 


r N ,fm 


A 




168.8 


172 


169 


1.255 


9.14 


3.74 


2.P1/2 


181.3 


185 


181 


1.350 


9.33 


3.72 


2si/ 2 


232 


239 


232 


1.745 


10.1 


3.64 


3^1/2 


254 


263 




1.920 


10.5 


3.60 



Footnote to the Table: the values of Z\ r correspond to the "naked nucleus" 
with the cut-off model (3'), Z cr are calculated with account of screening of the 
Coulomb field V(r) — —(/r by outward electrons (except of the if-shell, which 
is supposed to be ionized), £ cr = Z cr /137, and the values Z cr are calculated by 
the asymptotic formula (5). 

Note that the electric field E{r) near surface of a heavy nucleus, r s» rjy = 
1.2A 1 / 3 fm, is much larger than the "critical" or Schwinger field in QED [22]: 

E cr = m 2 e c 3 /eh = 1.32 • 10 16 V/cm, (6) 

and the Coulomb field of a heavy nucleus is equal to 

E{r) = Ze/r 2 , r > r N ; E{r N )/E cr = Za (l C /r N ) 2 - (7) 

So, E(r N ) w 1800 E cr for U nucleus (Z = 92, A = 238, r N = 7.45 fm) and 
E(r N ) « 

2200 E cr for Z = 172 (A = 2.6 Z, rjy = 9.2 fm). However, the static supercriti- 
cal field E(r) is strongly inhomogeneous and exists only in a small space region 
near rjy, therefore no e + e~ pairs can be produced by this field if Z < Z^. As is 
seen from eq.(7), the Coulomb field is larger than E cr only at distancesEl r < Iq. 

The same is true also for the field (17) of the massive nuclear density core, 
where 

E{r) « E max ■ r 2 - E max (—^\ , £»1, (8) 

the parameters £ and A are defined in Eq.(14) below and E{r) > E cr only near 
the core radius, 

r - Rc < \^E max /E cr ~ a-^Hc. (9) 

2 Note that E(l c ) = E cr if Z = a -1 = 137, and E(r) > E cr at r < \fZa l c . 
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Therefore the well-known formula [22] for pair production probability in homo- 
geneous electrostatic field, w oc (E/E cr ) 2 ■ exp(—irE cr /E), is not applicable in 
these cases. 

For Z > Z cr the vacuum becomes unstable with respect to production of 
e + e~-pairs. On account of the Pauli principle the number of produced pairs 
is determined by the number of discrete levels, which have descended into the 
lower continuum. Passing through the Coulomb barrier positrons go out to 
infinity, while electrons remain near the nucleus, partially screening its charge. 
Thus, a naked nucleus of supercritical charge Z > Z cr will envelop itself with 
an electron shell created out of the vacuum; we can call this shell "the vacuum 
shell" 

If Za ^> 1, the vacuum shell contains many electrons H and statistical ap- 
proach is necessary. The relativistic Thomas-Fermi equation [7,18,19] can be 
applied to calculate electron density n e {r). Let V(r) be the self-consistent po- 
tential for an electron, taking into account both the field of the nucleus and 
the average field created by other electrons of the vacuum shell. In WKB- 
approximation the electron momentum is 

p{r) = [{e-V{r)) 2 -m 2 e ] 1/2 , H = c=l (10) 

(in the WKB formula (10) the spin of electron is neglected, which is valid 
for large Z ^> 137). The vacuum shell of super- heavy nucleus is degenerated 
relativistic Fermi-gas with electron density 

n ^) = %T = ^(V 2 +2m e V) 3 / 2 , (11) 

where the value of P max follows from eq.(10) at e = — m e , since we are inter- 
ested only in the electrons that have dived into the continuum of the negative 
energy states. The spatial distribution of vacuum electrons is determined by 
the relativistic Thomas-Fermi equation 

Ay = -4vre 2 j JL (y 2 + 2m e Vf/ 2 - ra p (r) J (12) 

with the boundary conditions: V{oo) — (due to global charge neutrality of 
the system) and finiteness of V(0). Here n p (r) = n p 9(R C — r) is the proton 
density, n p = N p no/A^i 

0.25 m 3 , n — SA/^ttR 3 is the ordinary nuclear density, N p = Z is the number 
of protons and R c = Np^m^ 1 is the core radius. 

The density n e (r) of electrons is determined also by the Fermi energy con- 
dition on their Fermi momentum, Pf = P m ax'- 

El = [ (Pf ) + m 2 e Y' 2 - me - V(r) = 0, (13) 

which immediately leads to eq.(10). The equations for neutron, proton and 
electron densities have been integrated numerically [20]. 



3 The values of critical charge Z cr for highly excited atomic states were calculated in ref. [16]. 
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If N p e 3 >> 1, the electric field is concentrated in a narrow transition layer 
[7,18] of thickness ~ A s» 15 fm near r = R c » A therefore geometry reduces 
to the plane one. In the variables \ an d £ one nas □ 

V(r) = -{^n^'^x, £ = (r-R B )/\, 

X- 1 = 2(tt/3) 1/6 v / ^ n^ 3 sa ^ m w , (14) 

where A -C lc- 

Vlc-^k — nk- (14 ' } 
y/a ran 2b 

Therefore Eq.(12) becomes 

<Px/d? = X*-9{-t), x(~oo) = l, x(°o) = (15) 

and can be solved analitically [18]: 

f l-3[ 1 + 2- 1 / 2 sinh(a - eV3) £ < 0, 
X(0 = { (16) 
{ 2^(Z + b)-\ £>0, 

where 9(x] is the Heaviside step function and the integration constants are 

a = Arsh(llV2) = 3.439, b= -y/2 = 1.868. (16') 
Note that at £ < 0, i.e. inside the superheavy nucleus 

X (0 « 1 - 0.272 exp(€>/3) 1, U |» 1 (17) 
The electric field of the system 

/o5 \l/6 2 

^) = (^J v^x'co as) 

is damped exponentially inside the nucleus: E(£) oc exp(£\/3) as £ — > oo and 
^(0 ^ £ 2 — ^ in the outer region, r > R c . The field attains its maximal 
strength at the edge of super-charged nucleus 

TV?" ( 77? \ *" 

Emax = 0.95^ — « — #cr, f = i? c . (19) 

e V m e / 

So, E max s» 6000 E'er, which exceeds the characteristic field (6) in QED and is 
of the same order of magnitude as electric field (7) at surfaces of heavy nuclei 
with Za > 1. 



4 Note that the thickness of the transition layer A does not depend on values of the core 
radius R c and the mass number A, if R c 3> A. 
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In the region R c — r A the electric field is practically absent, and the 
electrically neutral plasma is formed inside the supercharged nucleus, where the 
densities n e and n p are equal and the potential is practically constant: 

9t^ 1/3 



V(r) « V(0) = - I — J Jn^ -1.92 m,r. (20) 

The uncompensated charge is situated in a layer of finite thickness ~ A near the 
edge of the nucleus, H/m^c <A< h/m e c. Though the formation of electrically 
neutral plasma inside a supercharged nucleus, Ze 3 3> 1, strongly diminishes the 
Coulomb energy of nucleus Q, but it remains positive and impedes the stability 
of such gigantic nuclei. So, the conclusion of ref. [18] is that nuclei with a mass 
number A ~ 10 4 - 10 6 are unstable due to the Coulomb repulsion of protons 
and can not exist in Nature. 

However, the situation may be changed considerably if one accounts a grav- 
itational attraction. Let us start with a simple qualitative estimate. The 
Coulomb energy is Ec ~ E^^R^X, which is mainly distributed within a thin 
shell of width A and radius R c 3> A. To ensure the stability of the system, the 
attractive gravitational energy of the shell (its mass m ~ MX/R C , M — Am n 
is mass of the core) 

£ _ GMm GM 2 \ _ Gm 2 n A 2 \ ^ 
R c R 2 R 2 

has to be larger than the repulsive Coulomb energy Ec- Since | £ gr |~ A 4 / 3 , 
while Ec ~ A 2 / 3 as A — > oo, a crossing | £ gr |= Ec necessarily exists: 

fl " Gm 2 n A 2 / 3 = (m n / mplanck ) 2 A 2 / 3 , (22) 



E c 



where mpi anc k — y/hc/G ~ 10 5 g is the Planck mass and m n ~ 10 24 g is the 
nucleon mass. So, | £ gr \> Ec at 

A > (m Planck /m n ) 3 » 10 57 , i? c = 1.2A 1 / 3 fm ~ 10 6 cm, (23) 

which are typical values for neutron stars. The Coulomb repulsion of protons, 
screened by relativistic electrons, is now balanced by gravitational forces. 

The more accurate derivation of gravitational and electrodynamical stability 
is based on the analytic solution (16) of the Thomas-Fermi equation. The 
Coulomb energy £ c and gravitational energy £ gr of the thin proton shell are 
[29] 

OO 

[E 2 R 3 (eV(0)) 3 f 

Ec = !^ dr= (3^)1/2 y[*(o]g = 



5 Due to the screening effect, the Coulomb energy = 3(N p e) 2 /5R c of a uniformly 

charged sphere (without screening) diminishes by 1.7 Ze 3 times of magnitude [16]. 
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0.345 



A 2jZ ( 



2/3 



V A 

GMm -3 Gmj 4/3 ( N p \ 1/3 m7I e 



(24) 



(25) 



where m is the mass of the layer and Gm^ = he (m n /mpi anc k) 2 . Hence, | E gr |> 
£ c at A > A R , 

^- 0,0390^) ( m ^ anCfc ) ' ( 26 ) 

which establishes a lower limit for the mass number A necessary for the stability 
of the massive nuclear density cores. 

However, besides the Coulomb energy £q, the kinetic energy of the degen- 
erated electronic Fermi-gas exists, 

Skin = n e (p)V = ^N p P Fl (27) 

where we took into account that the mean energy of particles in the degenerated 
relativistic Fermi-gas is (e) = (p) = jPf- The energy Skin also impedes the 
stability of the system and it should be compensated by the total gravitational 
energy of the core £q = —3GM 2 /5R C . A simple calculation shows that | Eg |> 
Ekin for 

N p \ ( mpia 



where c\ is a numerical constant of the order of unity. So, we again arrive 
at the condition, similar to Eq.(26), which is necessary for stability of massive 
nuclear density cores. Therefore it seems possible to formulate a consistent 
stable model of massive cores in terms of gravitational, strong, electromagnetic 
and weak interactions and quantum statistics. Certainly, many aspects of the 
problem of stability remain unsolved and further investigations are necessary. 

This work is based on the papers [18,29]. I would like to thank Profes- 
sor R.Rufhni, who initiated this work, and also V.Mur, L.Okun', M.Trusov, 
G.Vereshchagin, D.Voskresensky, and S.-S.Xue for valuable discussions and re- 
marks. The work was partially supported by ICRANet (Pescara, Italy) and 
Russian Foundation for Fundamental Research, project 07-02-01116. 
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